Machine learning detection of SARS-CoV-2 high-risk variants (preprint)

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has evolved many high-risk variants, resulting in repeated COVID-19 waves of pandemic during the past years. Therefore, accurate early-warning of high-risk variants is vital for epidemic prevention and control. Here we construct a machine learning model to predict high-risk variants of SARS-CoV-2 by LightGBM algorithm based on several important haplotype network features. As demonstrated on a series of different retrospective testing datasets, our model achieves accurate prediction of all variants of concern (VOC) and most variants of interest (AUC=0.96). Prediction based on the latest sequences shows that the newly emerging lineage BA.5 has the highest risk score and spreads rapidly to become a major epidemic lineage in multiple countries, suggesting that BA.5 bears great potential to be a VOC. In sum, our machine learning model is capable to early predict high-risk variants soon after their emergence, thus greatly improving public health preparedness against the evolving virus.

McAN: an ultrafast haplotype network construction algorithm (preprint)

Haplotype network is becoming popular due to its increasing use in analyzing genealogical relationships of closely related genomes. We newly proposed McAN, a minimum-cost arborescence based haplotype network con-struction algorithm, by considering mutation spectrum history (mutations in ancestry haplotype should be contained in descendant haplotype), node size (corresponding to sample count for a given node) and sampling time. McAN is two orders of magnitude faster than the state-of-the-art algorithms, making it suitable for analyzation of massive se-quences. Availability: Source code is written in C/C++ and available at https://github.com/Theory-Lun/McAN and https://ngdc.cncb.ac.cn/biocode/tools/BT007301 under the MIT license. The online web service of McAN is available at https://ngdc.cncb.ac.cn/ncov/online/tool/haplotype. SARS-CoV-2 dataset are available at https://ngdc.cncb.ac.cn/ncov/.

Global Cold-Chain Related SARS-CoV-2 Transmission (preprint)

The COVID-19 pandemic poses a great threat to human society. SARS-CoV-2 is mainly transmitted through social contact; however, it is highly debated whether cold-chain related transmission has occurred and can be identified in the epidemic areas of COVID-19. Here, we provide a new method and distinguish two transmission routes by detecting a lineage-specific reduction of SARS-CoV-2 mutation rate. After analyzing 1,610,125 SARS-CoV-2 genomic sequences, we find that two outbreaks in Xinfadi-Beijing and Auckland are cold-chain related and respectively caused by two mutation-dormant variants. A Dalian outbreak in July 2020 and a Yingkou outbreak ten months later are epidemiologically connected and derived from a cold-chain related variant. Mutation-dormant variants are detected during the spread of spike D614G variant and the Delta Variant of Concern. Cold-chain contaminations repeatedly caused by epidemiologically connected patients are also found and have resulted in infections. Moreover, the COVID-19 outbreak in Wuhan is likely to be cold-chain related. A systematic identification reveals that the frequency of cold-chain related transmission is in the order of magnitude of 0.1-10%. Our results indicate that that cold-chain related transmission is rare but happens globally.

Coronavirus GenBrowser for monitoring adaptive evolution and transmission of SARS-CoV-2 (preprint)

COVID-19 has widely spread across the world, and much research is being conducted on the causative virus SARS-CoV-2. To help control the infection, we developed the Coronavirus GenBrowser (CGB) to monitor the pandemic. CGB allows visualization and analysis of the latest viral genomic data. Distributed genome alignments and an evolutionary tree built on the existing subtree are implemented for easy and frequent updates. The tree-based data are compressed at a ratio of 2,760:1, enabling fast access and analysis of SARS-CoV-2 variants. CGB can effectively detect adaptive evolution of specific alleles, such as D614G of the spike protein, in their early stage of spreading. By lineage tracing, the most recent common ancestor, dated in early March 2020, of nine strains collected from six different regions in three continents was found to cause the outbreak in Xinfadi, Beijing, China in June 2020. CGB also revealed that the first COVID-19 outbreak in Washington State was caused by multiple introductions of SARS-CoV-2. To encourage data sharing, CGB credits the person who first discovers any SARS-CoV-2 variant. As CGB is developed with eight different languages, it allows the general public in many regions of the world to easily access pre-analyzed results of more than 132,000 SARS-CoV-2 genomes. CGB is an efficient platform to monitor adaptive evolution and transmission of SARS-CoV-2.

Structural classification of neutralizing antibodies against the SARS-CoV-2 spike receptor-binding domain suggests vaccine and therapeutic strategies (preprint)

The COVID-19 pandemic presents an urgent health crisis. Human neutralizing antibodies (hNAbs) that target the host ACE2 receptor-binding domain (RBD) of the SARS-CoV-2 spike1-5 show therapeutic promise and are being evaluated clincally6-8. To determine structural correlates of SARS-CoV-2 neutralization, we solved 8 new structures of distinct COVID-19 hNAbs5 in complex with SARS-CoV-2 spike trimer or RBD. Structural comparisons allowed classification into categories: (1) VH3-53 hNAbs with short CDRH3s that block ACE2 and bind only to "up" RBDs, (2) ACE2-blocking hNAbs that bind both "up" and "down" RBDs and can contact adjacent RBDs, (3) hNAbs that bind outside the ACE2 site and recognize "up" and "down" RBDs, and (4) Previously-described antibodies that do not block ACE2 and bind only "up" RBDs9. Class 2 comprised four hNAbs whose epitopes bridged RBDs, including a VH3-53 hNAb that used a long CDRH3 with a hydrophobic tip to bridge between adjacent "down" RBDs, thereby locking spike into a closed conformation. Epitope/paratope mapping revealed few interactions with host-derived N-glycans and minor contributions of antibody somatic hypermutations to epitope contacts. Affinity measurements and mapping of naturally-occurring and in vitro-selected spike mutants in 3D provided insight into the potential for SARS-CoV-2 escape from antibodies elicited during infection or delivered therapeutically. These classifications and structural analyses provide rules for assigning current and future human RBD-targeting antibodies into classes, evaluating avidity effects, suggesting combinations for clinical use, and providing insight into immune responses against SARS-CoV-2.

The global landscape of SARS-CoV-2 genomes, variants, and haplotypes in 2019nCoVR (preprint)

On 22 January 2020, the National Genomics Data Center (NGDC), part of the China National Center for Bioinformation (CNCB), created the 2019 Novel Coronavirus Resource (2019nCoVR), an open-access SARS-CoV-2 information resource. 2019nCoVR features a comprehensive integration of sequence and clinical information for all publicly available SARS-CoV-2 isolates, which are manually curated with value-added annotations and quality evaluated by our in-house automated pipeline. Of particular note, 2019nCoVR performs systematic analyses to generate a dynamic landscape of SARS-CoV-2 genomic variations at a global scale. It provides all identified variants and detailed statistics for each virus isolate, and congregates the quality score, functional annotation, and population frequency for each variant. It also generates visualization of the spatiotemporal change for each variant and yields historical viral haplotype network maps for the course of the outbreak from all complete and high-quality genomes. Moreover, 2019nCoVR provides a full collection of SARS-CoV-2 relevant literature on COVID-19 (Coronavirus Disease 2019), including published papers from PubMed as well as preprints from services such as bioRxiv and medRxiv through Europe PMC. Furthermore, by linking with relevant databases in CNCB-NGDC, 2019nCoVR offers data submission services for raw sequence reads and assembled genomes, and data sharing with National Center for Biotechnology Information. Collectively, all SARS-CoV-2 genome sequences, variants, haplotypes and literature are updated daily to provide timely information, making 2019nCoVR a valuable resource for the global research community. 2019nCoVR is accessible at https://bigd.big.ac.cn/ncov/.

SARS-CoV-2 infects human pluripotent stem cell-derived cardiomyocytes, impairing electrical and mechanical function (preprint)

Global health has been threatened by the COVID-19 pandemic, caused by the novel severe acute respiratory syndrome coronavirus (SARS-CoV-2)1. Although considered primarily a respiratory infection, many COVID-19 patients also suffer severe cardiovascular disease2-4. Improving patient care critically relies on understanding if cardiovascular pathology is caused directly by viral infection of cardiac cells or indirectly via systemic inflammation and/or coagulation abnormalities3,5-9. Here we examine the cardiac tropism of SARS-CoV-2 using human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs) and three-dimensional engineered heart tissues (3D-EHTs). We observe that hPSC-CMs express the viral receptor ACE2 and other viral processing factors, and that SARS-CoV-2 readily infects and replicates within hPSC-CMs, resulting in rapid cell death. Moreover, infected hPSC-CMs show a progressive impairment in both electrophysiological and contractile properties. Thus, COVID-19-related cardiac symptoms likely result from a direct cardiotoxic effect of SARS-CoV-2. Long-term cardiac complications might be possible sequelae in patients who recover from this illness.

A simplified cell-based assay to identify coronavirus 3CL protease inhibitors (preprint)

We describe a mammalian cell-based assay capable of identifying coronavirus 3CL protease (3CLpro) inhibitors without requiring the use of live virus. By enabling the facile testing of compounds across a range of coronavirus 3CLpro enzymes, including the one from SARS-CoV-2, we are able to quickly identify compounds with broad or narrow spectra of activity. We further demonstrate the utility of our approach by performing a curated compound screen along with structure-activity profiling of a series of small molecules to identify compounds with antiviral activity. Throughout these studies, we observed concordance between data emerging from this assay and from live virus assays. By democratizing the testing of 3CL inhibitors to enable screening in the majority of laboratories rather than the few with extensive biosafety infrastructure, we hope to expedite the search for coronavirus 3CL protease inhibitors, to address the current epidemic and future ones that will inevitably arise.

Single-cell analysis of human trophoblast stem cell specification reveals activation of fetal cytotrophoblast expression programs including coronavirus associated host factors and human endogenous retroviruses (preprint)

The human placenta is increasingly a focus of research related to early child development and the impact of maternal hyperimmune states. The ability to model human trophectoderm disease states from human pluripotent stem cells, the nature of human pluripotent stem cell potency and the mechanisms regulating human trophectoderm specification remains poorly understood. Recent work suggests that only the naive state can give rise to trophectoderm and that primed iPSC generate mixed amnionic and mesoderm lineages. Here we identify conditions that efficiently drive the specification of primed iPSC to trophectoderm, named Trophoblast Stem Cell (TSC). iPS-derived-TSC share transcriptional, morphological and functional characteristics with human in vivo cytotrophoblasts including activation of human endogenous retroviruses, expression of COVID-19 associated host factors and generation of multinucleated syncytiotrophoblasts with a large fusion index. At high densities in 5% O2, iPS-derived-TSC form villi-like structures and express extravillous and syncytiotrophoblast proteins HCG-{beta} and HLA-G. Using temporal single cell RNAseq, we define the molecular changes associated with specification under three separate conditions: 1) BMP4, 2) BMP4 and inhibition of WNT, 3) activation of EGF and WNT, inhibition of TGFbeta, HDAC and ROCK signaling (named TSC). With 9,821 high-quality single cell transcriptomes, we find that BMP4 gives rise to mesenchymal cells while TS conditions lacking exogenous BMP4 generate a stable proliferating cell type that is highly similar to six week placenta cytotrophoblasts. TFAP2A programs the specification of primed iPS cells to TSC without transitioning through a naive state. TSC specification independent of exogenous BMP4 will allow for robust and reproducible studies of the cytotrophoblast component of human placenta.